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We suggest the idea, supported by concrete calculations within chiral models, that the critical 
endpoint of the phase diagram of Quantum Chromodynamics with three colors can be detected, 
by means of Lattice simulations of grand-canonical ensembles with a chiral chemical potential, /is, 
conjugated to chiral charge density. In fact, we show that a continuation of the critical endpoint of 
the phase diagram of Quantum Chromodynamics at finite chemical potential, /i, to a critical end 
point in the temperature-chiral chemical potential plane, is possible. This study paves the way of 
the mapping of the phases of Quantum Chromodynamics at finite fi, by means of the phases of a 
fictitious theory in which y, is replaced by /is. 
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I. INTRODUCTION 



The critical endpoint, CP, of Quantum Chromodynam- 
ics (QCD) [l[ is one of the most important aspects of 
the phase diagram of strongly interacting matter. It 
is thus not surprising that an intense experimental ac- 
tivity is nowadays dedicated to the detection of such a 
point, which involves the large facilities at RHIC and 
LHC; moreover, further experiments are expected after 
the development of FAIR at GSI. Several theoretical sig- 
natures of CP have been suggested 0, 0]. Despite the 
importance of CP, a firm theoretical evidence of its exis- 
tence is still missing. In fact, the sign problem makes the 
Lattice Monte Carlo simulations difficult, if not impossi- 
ble, in the large baryon-chemical potential (/i) region for 
N c = 3 ||, see for a recent review. Therefore, it has 
not yet been possible to prove unambiguously the exis- 
tence and the location of CP starting from first principles 
simulations of grand-canonical ensembles. The strong 
coupling expansion of Lattice QCD [f| 0] seems promis- 
ing. Even more, the predictions of effective models are 
spread in the T — /i plane, see for example 

An interesting overcoming of the sign problem for the 
quest of CP is offered by analytic continuation of data ob- 
tained at imaginary chemical potential, \ii (lfj|-[l^|. Re- 
cent promising analysis shows that it might be possible 
to continue the critical line from the region of imaginary 
/i to that of real /i (l3j |. pinning down the critical point. 
Another fruitful approach is given by simulations at fi- 
nite isospin chemical potential, S/i, see for example [lil - 
0]|. Even in this critical endpoint there appears 

before the transition to the pion condensed phase. The 
latter point has been overlooked and it has not yet been 
detected by mean field model calculations; hence it cer- 
tainly deserves further study. It is also worth to cite 
the possibility to perform simulations in canonical, rather 
then grand-canonical, ensembles. Preliminary results in 
this direction have been presented recently in [l7|. On 
the purely theoretical side, it has been suggested very 



recently [l8| that the use of orbifold equivalence in the 
large N c approximation of QCD can lead to relations be- 
tween the coordinates of CP at finite chemical potential, 
with those at finite isospin chemical potential. 

In this Article, we suggest a new, theoretical way to de- 
tect the CP, by means of Lattice simulations with N c = 3, 
which can be considered as an alternative to the \xi tech- 
nique. In order to accomplish this important program, we 
suggest to simulate QCD with a chiral chemical potential, 
/i5, co njug ated to the chiral charge density, = n^-ni, 
see [19-22] for previous studies. Our idea, supported by 
concrete calculations within microscopic effective mod- 
els, is that CP can be continued to a critical endpoint 
at fi 5 ^ and fi = 0, that we denote by CP5, the lat- 
ter being accessible to iV c = 3 Lattice QCD simulations 
of grand-canonical ensembles [n|. Therefore, the detec- 
tion of the former endpoint via Lattice simulations, can 
be considered as a signal of the existence of the latter. 
To facilitate exposition, we introduce the symbol W5 to 
denote the world with /i = 0, [15 ^ 0. On the other 
hand, we will use the symbol W to denote the world with 
/i5 = 0, and which corresponds to the physical universe. 

The model calculations, in particular the ones based 
on the Nambu-Jona-Lasinio model with the Polyakov 
loop [23| (PNJL model in the following) with tree level 
coupling among chiral condensate and Polyakov loop [24j , 
give numerical relations among the coordinates of CP5 
and those of CP. In particular, the critical temperature 
turns out to be almost unaffected by the process of con- 
tinuation; the critical value of the chemical potential, [i c , 
on the other hand turns out to be almost half of the 
critical chiral chemical potential, fj,§ c . 

Before discussing our results, it is important to spend 
some word more about the chiral chemical potential. In 
particular, we are aware that W5 should be considered 
as a fictional universe. As a matter of fact, ^ cannot be 
considered as a true chemical potential because, in the 
confinement phase, the chiral condensate (qq) mixes left- 
and right-handed components of the quark field, leading 
to non-conservation of n$. This statement is true also 
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in the Quark-Gluon-Plasma phase, where the chiral con- 
densate is much smaller than its value in the confinement 
phase; in this case, the non-conservation of 72,5 is much 
softer, and mainly due to the bare quark mass, m <C T 
with T corresponding to the temperature of the heath 
bath in which the fields live. Therefore, the point of view 
that we adopt in this Article is to consider /15 as a mere 
mathematical artifice. However, W5 with N c = 3 can be 
simulated on the Lattice. Therefore, for the continuity 
property cited above, it is worth to study it by grand- 
canonical ensemble simulations: it might furnish an evi- 
dence of the existence of the critical endpoint in the real 
world. Furthermore, once CP5 is detected, it might be 
possible to make use of Lattice simulations to detect in- 
homogeneous phases [2M [26{ which could develop around 
CP5, as a continuation of those which develop at CP. The 
results would then be of vital importance to understand, 
for example, the inner structure of compact stellar ob- 
jects. For the aforementioned reasons, this study is very 
far from being of purely academic or theoretical interest. 



II. CHIRAL MODELS 

Because of its non-perturbative nature, we cannot 
make first principles calculations within QCD in the 
regimes to which we are interested in, namely moder- 
ate T, fi and fj.5 . Hence we need to rely on some effective 
model, which is built in order to respect (at least some of) 
the symmetries of the QCD action. To this end, we make 
use of the celebrated Quark- Meson model [27j , and of the 
Nambu- Jona-Lasinio model f2S || (s ee [29j | for reviews) im- 
proved with the Polyakov loop [23], dubbed PNJL model, 
which have been used many times in recent years to de- 
scribe successfully the thermodynamics of QCD with two 
and two-plus-one flavors, see [24 . l3CH37l | and references 
therein. They are interesting because they allow for a 
self-consistent description of spontaneous chiral symme- 
try breaking; even more, in the case the model is im- 
proved with the Polyakov loop, it allows for a simulta- 
neous computation of quantities sensible to confinement 
and chiral symmetry breaking. 

In this Section we describe both of the models. We 
firstly discuss the Quark- Meson model that we use in our 
computation, which is simpler to implement since take its 
simplest version, without the complications due to the fi- 
nite value of the quark masses and to the Polyakov loop. 
This is an useful preparation to the more important case 
of the PNJL model, which is more trustable quantita- 
tively since it is tuned to reproduce Lattice data at zero 
and imaginary chemical potential. 



A. Quark-Meson Model 

The Quark-Meson (QM) model consists of the 0(4) 
linear sigma model coupled to dynamical quarks via a 
Yukawa-type interaction. The lagrangian density is given 



by 

+ \{d,a) 2 + l -(d^f-U(a^) . (1) 

In the above equation, q corresponds to a quark field 
in the fundamental representation of color group SU (3) 
and flavor group SU{2); besides, a and 7r correspond to 
the scalar singlet and the pseudo-scalar iso-triplet fields, 
respectively. We have a introduced chemical potential 
for the quark number density, fj,, and a pseudo-chemical 
potential conjugated to chirality imbalance, /U5. 

We explain in some detail the physical meaning of the 
latter. The quantity conjugated to u 5 , namely the chiral 
charge density, is given by 115 = ur — Ul, and represents 
the difference in densities of the right- and left-handed 
quarks. At finite /U5, a chirality imbalance is created, 
namely TI5 ^ 0. For example, in the massless limit and 
at zero baryon chemical potential one has [l9| 



A 

3tt 2 



(2) 



If quark mass (bare or constituent) is taken into account, 
the relation 1x5(115) cannot be found analytically in the 
general case, and a numerical investigation is needed, see 
for example [2(3]. The imbalance of chiral density can 
be created by instanton/sphaleron transition in QCD, 
see and references therein. As a matter of fact, 
non perturbative background gluon configurations with 
nonzero winding number, Qw, change the chirality of 
quarks according to the Ward identity: 



(in 5 
~~dt~ 



9 2 N f 
16tt 2 



d 3 xF^Pr 



Qw 



(3) 



As a consequence, the addition of /15 to the la- 
grangian density of the chiral models mimics the instan- 
ton/sphaleron induced chirality transitions. 

The potential U describes tree-level interactions among 
the meson fields. In this Article, we restrict ourselves to 
the QM model in the chiral limit for simplicity, and take 



(4) 



which is invariant under the chiral group. 

We work in the one-loop approximation, which 
amounts to consider mesons as classical fields, and in- 
tegrate only over fermions in the generating functional 
of the theory to obtain the Quantum Effective Potential 
(QEP). In the integration process, the meson fields are 
fixed to their classical expectation values, (tt) = and 
(a) 7^ 0. The physical value of (a) will be then deter- 
mined by minimization of the QEP. The field a has the 
quantum numbers of the QCD chiral condensate, (qq). 
Hence its non-vanishing expectation value breaks chiral 
symmetry spontaneously, mimicking the chiral symmetry 
breaking of the QCD vacuum. 
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The one-loop QEP in presence of /Z5 has been discussed 
several times [19l - l21| . The addition of the \i— dependence 
is a textbook matter. The final result is 

d 3 p 



N c N f 



d 3 p 
d 3 p 



where 



(Ip|s-M5) 5 



, (5) 



(6) 



corresponds to the pole of the quark propagator, and 
m q = ga is the constituent quark mass; finally, the index 
s denotes the helicity projection. 

In right hand side of the first line of Equation (|5j) the 
momentum integral corresponds to the vacuum quark 
fluctuations contribution to the QEP. It is divergent, and 
it gives a contribution at T = and /Z5 ^ 0. For our 
purposes, it is enough to treat the model as an effec- 
tive description of the infrared regime of QCD. There- 
fore we treat the divergence phenomenologically, intro- 
ducing a momentum cutoff, M, in the vacuum term. 
This is equivalent to introduce a momentum-dependent 
quark mass, which is nonzero and constant for momenta 
lower than the cutoff, and zero for larger values of mo- 
menta, thus realizing a rough approximation of the effec- 
tive (momentum-dependent and ultraviolet suppressed) 
quark mass of full QCD [3a |. 

The parameters of the model are tuned as in [9[ . They 
are chosen in order to satisfy the requirements that a — 
f„ in the vacuum, 



and m„ 



dV 

700 MeV, where 







d 2 V 



da 2 



(7) 



(8) 



Moreover, the constituent quark mass in the vacuum is 
fixed to the value m q = 335 MeV, which allows to fix 
the numerical value of g = m q /f„ with / w = 92.4 MeV. 
Finally, the ultraviolet cutoff is taken M — 600 MeV. 
This procedure fixes A = 2.73 and v 2 = -(617.7 MeV) 2 . 



B. PNJL model 

In the PNJL model, quark propagation in the medium 
is described by the following lagrangian density: 



£ = q(ij> 1 D li -m,)q + £ I ; 



(9) 



here q is the quark Dirac spinor in the fundamental repre- 
sentation of the flavor SU(2) and the color group; r cor- 
respond to the Pauli matrices in flavor space. A sum over 
color and flavor is understood. The covariant derivative 
embeds the QCD coupling with the background gluon 
field which is related to the Polyakov loop, see below. 
Furthermore, we have defined 



Ci = G (qq) + (iq^rqf 



(10) 



One advantage to use the PNJL model is that it has 
access to the expectation value of the Polyakov loop, that 
we denote by L, which is sensible to confinement or de- 
confinement properties of a given phase. In order to com- 
pute L we introduce a static, homogeneous and Euclidean 
background temporal gluon field, Aq = iA^ — i\ a A% : 
coupled minimally to the quarks via the QCD covariant 
derivative [23j. Then 



iTr c exp {i/3\ a Al) 



(11) 



where /3 — l/T. In the Polyakov gauge, which is con- 
venient for this study, Aq — + i^9,4> s ] moreover, for 
simplicity we take L = from the beginning as in [30j, 
which implies Af = 0. We have then 



1 + 2 cos(/30) 



(12) 



For our purpose, that is mainly the location of the critical 
endpoint, we expect that the approximation L = is 
sufficient; as a matter of fact, we have verified that in this 
case we reproduce, at finite /z, the location of the critical 
endpoint obtained in (24[, where the same parametriza- 
tion of the model is used, and where L 7^ L' from the 
beginning. 

In our computation we follow the idea implemented 
in [24], which brings to a Polyakov- loop-dependent cou- 
pling constant: 



G = g [1 - ai Ltf - a 2 (L 3 + (L*) 3 )] 



(13) 



The ansatz in the above equation was inspired by [39|, |4fJ 
in which it was shown explicitly that the NJL vertex can 
be derived in the infrared limit of QCD, it has a non- 
local structure, and it acquires a non-trivial dependence 
on the phase of the Polyakov loop. We refer to [24| for 
a more detailed discussion. This idea has been analyzed 
recently in (4l| . where the effect of the confinement or- 
der parameter on the four-fermion interactions and their 
renormalization-group fixed-point structure has been in- 
vestigated. The numerical values of oti and a-i have been 
fixed in [24[ by a best fit of the available Lattice data 
at zero and imaginary chemical potential of Ref. [42|, |43[ . 
In particular, the fitted data are the critical temperature 
at zero chemical potential, and the dependence of the 
Roberge- Weiss endpoint on the bare quark mass. The 
best fit procedure leads to a.\ — a 2 = a = 0.2 ± 0.05, 
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within the hard cutoff regularization scheme, which is 
the same scheme that we adopt in this Article. 

In the one-loop approximation, the effective potential 
of this model is given by 



d 3 p 



V=U{L 1 L\T)+ U --N C N S E / 7^ 

s=±l J V 



Mi 



(14) 



where 



F_ = 1 + SLe-^'-ri + 3Lh- 2 ^'-^ + e'^"-^ , 

(15) 

F+ = 1 + 3L t e-' 9( " s+ ' i) + 3Le- mu » + ^ + e -80(".+/O , 

(16) 

denote the statistical confining thermal contributions to 
the effective potential; u> 3 is still given by Equation ([B]), 
with m q = m — 2a. Once again the vacuum fluctua- 
tion term is regularized by means of a ultraviolet cutoff, 
that we denote by M. The relation between the chiral 
condensate and a in the PNJL model is a = 2G(qq). 

We notice that in this case we take quarks with a finite 
bare mass, which will be fixed by requiring that the pion 
mass in the vacuum is in agreement with its experimen- 
tal value. We also notice that the PNJL model consid- 
ered here, which is dubbed Extcnded-PNJL in [24|, has 
been tuned in order to reproduce quantitatively the Lat- 
tice QCD thermodynamics at zero and imaginary quark 
chemical potential. Hence, it represents a faithful de- 
scription of QCD, in terms of collective degrees of free- 
dom related to chiral symmetry breaking and deconfine- 
ment. 

The potential term U in Eq. (|14p is built by hand in 
order to reproduce the pure gluonic lattice data with 
N c = 3 [30| . We adopt the following logarithmic form, 



U{L,L,T] = T i l [ -^±lL 

+ b{T) ln[l - 6LL + 4(Z 3 + L 3 ) - 3(LL) 2 ] 



(17) 



with three model parameters (one of four is constrained 
by the Stefan-Boltzmann limit), 
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FIG. 1. (Color online). Upper panel: Phase diagram of the 
QM model in the ^5 — T plane. The dot-dashed line corre- 
sponds to the second order chiral phase transition; the solid 
line denotes the first order phase transition. The orange dot is 
the critical endpoint. The scale T° = 174.1 MeV corresponds 
to the critical temperature at ^5 = 0. Lower panel: Phase 
diagram of the PNJL model. Lines convention is the same as 
in the upper panel. The scale T° = 173.9 MeV corresponds 
to the critical temperature at fj,s = 0. 



gauge theory. In absence of dynamical fermions one has 
To = 270 MeV. However, dynamical fermions induce a 
dependence of this parameter on the number of active 
flavors [34|. For the case of two light flavors to which 
we are interested here, we take To = 190 MeV as in fzij . 
Also for the remaining parameters we follow [24j and take 
M = 631.5 MeV, m = 5.5 MeV and G = 5.498 x 10~ 6 
MeV~ 2 . 



,I(T) «n f,/| ( + fl 2 Cj7 



(18) 



The standard choice of the parameters reads ao = 3.51, 
ai = —2.47, a2 = 15.2 and 63 = —1.75. The parameter 
To in Eq. (|17p sets the deconfinement scale in the pure 



III. CRITICAL ENDPOINT AT ZERO 
CHEMICAL POTENTIAL 

In Figure [TJ we plot the phase diagram of the chiral 
models in the /15 — T plane, for the case fj, = 0. It is 
obtained by a minimization procedure of the full poten- 
tial ([5]) for the QM model, and (fL4")l for the case of the 
PNJL model. 
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In the case of the QM model, the critical tempera- 
ture T c is identified with the one at which (a)x=T c = 0. 
For the PNJL model, since chiral symmetry is broken 
explicitly by the quark mass and the phase transitions 
are replaced by crossovers, we identify the critical tem- 
perature with that at which dL/dT is maximum. We 
have checked that the latter deviates from that at which 
\da / dT\ is maximum only of a few MeV, in the whole 
range of parameters studied. Even in this case, with an 
abuse of nomenclature, we dub the pseudo-critical lines 
as second order and first order, as in the case of the QM 
model. It is clear from the context that, whenever we 
talk about the PNJL model, the term second order tran- 
sition has to be taken as a synonym of smooth crossover; 
similarly, the term first order transition is a synonym of 
discontinuous jump of the order parameters. 

In the case of the QM model, the data about the con- 
densates as a function of temperature atjk ^ have 
been presented already in the literature [21], therefore 
their reproduction in this Article would be just of aca- 
demic interest; hence we skip this step here. Our focus is 
the discussion of the critical endpoint and of its evolution 
from W5 to W; hence we focus on CP5 and CP. 

In both panels of Figure [TJ the grey dashed line cor- 
responds to the second order chiral phase transition in 
the case of the QM model, or to a smooth crossover in 
the case of the PNJL model. The solid line, on the other 
hand, denotes the first order transition. The dot corre- 
sponds to CP5. In the following, we label the coordinates 
of CP5 by (n5c,T5 C ). For the case of the QM model, 
since we do not have any information about confining- 
deconfining property of a given phase, we can label the 
phases of the model only in terms of the chiral symme- 
try. We call the phase below the critical line as the chiral 
symmetry broken phase; similarly, above the critical line, 
chiral symmetry is restored, hence we call this phase as 
the chiral symmetry restored phase. It is then natural 
that CP5 in this case is a chiral critical endpoint. 

For the case of the PNJL model, on the other hand, we 
have access to the chiral condensate and to the Polyakov 
loop expectation value. As a consequence, we can label 
the phases of the model in terms both of confining prop- 
erties, and of chiral symmetry. In the model at hand, 
because of the entanglement in Equation (p~3|) . the de- 
confinement and chiral symmetry restoration crossovers 
take place simultaneously both at value of fj,§. This is 
proved by our data about a and L, see Figure [5] in which 
we plot the chiral condensate (upper panel) and the ex- 
pectation value of the Polyakov loop (lower panel) as a 
function of temperature, for several values of /15. In (l9j 
the two crossovers joined only for /15 > Our results 
are a natural consequence of the entan glem ent vertex in 
Equation (fl~5)) , which was neglected in [2(| • 

It is useful to stress that in the PNJL model we are 
discussing crossovers. Due to the crossover nature of 
the phenomena, a unique definition of the critical tem- 
peratures is not available; in this article we have used 
the definition which is easiest to implement, namely the 
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FIG. 2. (Color online) . Upper panel: Chiral condensate in the 
PNJL model, normalized to the value at zero temperature, as 
a function of temperature for several values of ^5. Black solid 
line corresponds to the case ^5 = 0; orange dashed line to 
lis = 200 MeV; grey dotted line to fi 5 = 300 MeV; indigo 
dot-dashed line to ^5 = 400 MeV. Lower panel: Expectation 
value of the Polyakov loop as a function of temperature, for 
several values of /15 . Dashing and color convention is the same 
as in the upper panel. 



study of the peaks of the derivatives of the order param- 
eters. We dub these quantities as effective susceptibili- 
ties. Within numerical error, we find that the peaks for 
the effective susceptibilities of the Polyakov loop and of 
the chiral condensate coincide in temperature, within few 
MeV. We cannot exclude that using different definitions 
for the pseudo-critical temperatures, like the identifica- 
tion of the crossovers with the peaks of true suscepti- 
bilities, the former can be shifted of some MeV, leading 
eventually to a larger split of the deconfinement and the 
crossover temperatures. However, the qualitative picture 
should not be modified drastically, as previous studies 
within the PNJL model have shown. 

Because of this peculiarity of the PNJL model, at 
the pseudo-critical line both deconfinement and chiral 
restoration crossovers take place. Hence the region below 
the pseudo-critical line is characterized by confinement 
and spontaneous breaking of chiral symmetry; we label 
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this phase as the confinement phase. On the other hand, 
the phase above the critical line is identified with the 
Quark-Gluon-Plasma phase. In this case, CP5 is both 
chiral and deconfinement critical endpoint. 

For what concerns the coordinates of CP5, in the case 
of the QM model we find 

f^, ^ - (2.16,0.78) , CP 5 (QM) , (19) 

where Xj? = 174.1 MeV is the chiral symmetry restoration 
temperature at /J, = //g = 0. Moreover, for the PNJL 
model we find 



M5c T c 
TP'TP 



(1.73,0.96), CPs (PNJL) , (20) 



where T° = 173.9 MeV is the deconfinement temperature 
at fi = /is = 0. 



IV. CRITICAL ENDPOINT AT FINITE 
CHEMICAL POTENTIAL 
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Next we turn to discuss the more general case with 
both /is and /i different from zero. Our scope is to show 
that, at least within the models, CP naturally evolves 
into CP5. Hence the latter, if detected on the Lattice, can 
be considered as the benchmark of the former. In partic- 
ular the PNJL model, which is in quantitative agreement 
with the Lattice at zero chemical potential, gives a nu- 
merical relation among the coordinates of CP and CP5, 
which might be taken as a guide to estimate the coordi- 
nates of CP in QCD, once CP5 is detected. 

In Figure [3] we collect our data on the critical point of 
the phase diagram in the /i — /X5 — T space, in the case 
of the PNJL model (for the QM model we obtain similar 
results) . The orange solid line is the union of the critical 
points computed self-consistently at several values of fi: 
at any value of fx, a point on the line corresponds to the 
critical point of the phase diagram in the /is — T plane. 
Thus the line pictorially describes the evolution of the 
critical point of the chiral model at hand, from CP to 
CP5. In Figure |4]we plot a projection of Fig. [3]onto the 
fi — /is plane, for the PNJL model. The indigo solid line 
corresponds to the /is-coordinate of the critical endpoint. 
The critical temperature is not so much affected when we 
continue CP5 to CP (we measure a change approximately 
equal to the 3%), therefore the projection in the \i — T 
plane is redundant. 

At /is = 0, the critical point is found at the following 
coordinates: 




= (0.92,0.93) , CP (PNJL) , (21) 



which correspond to /i c « 160 MeV and T c « 165 MeV, 
in agreement with the results of [24| , Similarly for the 
case of the QM model we find 




= (1.58,0.49), CP (QM) . (22) 



FIG. 3. (Color online). Evolution of the critical endpoint in 
the fi — /is — T space, for the PNJL model. 




50 100 150 

fi (MeV) 

FIG. 4. (Color online). Projection of Fig. [3] onto the /1 — ^5 
plane, for the PNJL model. The solid line corresponds to the 
/^-coordinate of the critical endpoint. 



The natural question which arises is: suppose a grand- 
canonical ensemble simulation finds CP5; then, how can 
the coordinates of the latter point help to locate CP? We 
can answer to this important question within the PNJL 
model. Our numerical results, Equations (|2T)]) and (j2"T|) . 
suggest the following relations: 

— « 0.53 , « 0.97 , (PNJL) . (23) 

A*5c 2 5c 

Within the QM model we obtain similar relations, but 
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in our opinion, those in (I23[) are more trustable quanti- 
tatively because the PNJL model has been tuned to be 
in quantitative agreement with Lattice data at zero as 
well as imaginary chemical potential [2~H ]. a characteris- 
tic which is not satisfied by the QM model used here. 

The model predictions (|23|) relate the coordinates of 
CP to those of CP5. In particular, it is interesting that 
the critical temperature is almost unchanged in the con- 
tinuation of CP to CP5. Of course, since these results are 
deduced by a model, it is extremely interesting and im- 
portant to study how Equation (|23p is affected by varying 
parameters like the bare quark mass, or the number of 
active flavors. This observation opens the possibility to 
develop further non-academic model studies of the prob- 
lem that we discuss in this Article. 



V. DISCUSSION 

In this Section we summarize briefly the results ob- 
tained, and discuss their conceptual relevance, and po- 
tential applications to the Lattice as well. 

Our main goal is to show that the continuation of the 
critical endpoint of the QCD phase diagram, CP, to a 
fictitious critical endpoint, CP5, belonging to a phase 
diagram in the /15 — T plane, is reasonable. Here /U5 cor- 
responds to the chiral chemical potential, which is conju- 
gated to the chiral density imbalance, n.5 = ur — n^. As 
we have already stressed, strictly speaking /j, 5 should be 
regarded as a pseudo-chemical potential: quark conden- 
sate in the confinement phase induces a mix among left- 
and right-handed quark field components. As a conse- 
quence, the conjugated quantity to /X5, namely the chiral 
charge density, is not conserved. Hence, it is legitimate 
to consider ^5 as a mere mathematical artifact, and the 
world in which [i is replaced by that we have bap- 
tized W5 , as an artificial universe, not necessarily related 
to the physical world. This is the point of view that we 
adopt in this Article. 

However, even accepting this minimalist point of view, 
W5 has the quality that it can be simulated on a Lat- 
tice. Indeed, the sign problem which affects simulations 
of three color QCD at finite chemical potential, does not 
affect QCD in W5 19]. Thus, it is possible to check 
wether CP5 there exists or not, within first principles 
calculations. If Lattice simulations find CP5, our result 
suggests the simplest interpretation: CP5 would be noth- 
ing but the continuation of the critical point which there 
exists in QCD. Hence we suggest to interpret CP5, if 
found, as a signal of the existence of the critical point in 
QCD. The coordinates of CP are related to those of CP5 
by Equation (12511 . 

Furthermore, it is extremely intriguing, theoretically 
speaking, that the critical endpoint CP5 that we find in 
this study, which confirms the findings of previous studies 
obtained within different models, is the analytic continu- 
ation of the critical endpoint, CP, of the theory at finite 
quark chemical potential. 



As anticipated, previous studies [2CJ, |21| have already 
discussed the possible appearance of CP5. In particular, 
in 20] a PNJL model in a strong magnetic background, 
without entanglement vertex [see Equation (|13[) ] has been 
considered. In that context, the chiral chemical potential 
was introduced to mimic the presence of chirality imbal- 
ance induced by instantons and sphalerons transitions in 
the hot QCD medium. Quantitatively, our results on the 
critical endpoint are in agreement with those of [20| . 

In [2l|, the phase diagram in the ^5 — T plane has 
been computed within the QM model improved with the 
Polyakov loop. The models considered here and in [2l[ 
are different: the zero point energy is not taken into ac- 
count in 2l|, and on the other hand we do not include 
the Polyakov loop in our QM model calculations. For 
these reasons, it is not necessary to have a quantitative 
agreement among our results and those of [21| . On the 
contrary, it is important to notice that the qualitative 
picture is the same, namely the existence of a critical 
endpoint in the ^5 — T plane. 

In both of the aforementioned studies, the quark chem- 
ical potential fi has not been taken into account. There- 
fore, the possibility to continue CP5 to CP, which is the 
main idea of the study illustrated in this Article, was not 
considered in those references. 



VI. CONCLUSIONS AND OUTLOOK 

In this Article, we have suggested the possibility of con- 
tinuation of the critical endpoint of the phase diagram of 
N c = 3 QCD, CP, to a critical endpoint dubbed CP5 be- 
longing to a fictional world, W5, in which the quark num- 
ber chemical potential is replaced by a pseudo-chemical 
potential, /Z5 conjugated to the chiral charge density, n.5. 
The universe VV5 has the merit that it can be simulated 
on the Lattice [19| for N c = 3. This suggestion is based 
on concrete calculations within chiral models. In partic- 
ular, we have used the PNJL model with entanglement 
vertex, introduced in (24|, which offers a description of 
the QCD thermodynamics in terms of collective degrees 
of freedom, which is in quantitative agreement with Lat- 
tice data at zero and imaginary chemical potential. 

Our main idea is that simulations in W5 might reveal 
the existence of a critical endpoint, CP5, in the phase 
diagram. Then, this critical point might be interpreted as 
the continuation of the critical point which is expected to 
belong to the phase diagram of real QCD, because of the 
continuity summarized in Fig. [5J Hence it would be an 
indirect evidence of the existence of the critical point in 
real QCD. Moreover, numerical predictions of the model, 
which connect CP to CP5, are in Equation (|23p . For 
these reasons, the interest of the present study is very 
far from being only academic or purely theoretical. Our 
result paves the way of the mapping of the phases of 
Quantum Chromodynamics at finite [A, by virtue of the 
phases of a fictitious theory in which /i is replaced by ^5 . 

In our calculations there are some factors that we have 
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not included for simplicity, and that might be interesting 
to include in more complete calculations (massive quarks, 
vector interactions, just to cite a couple of examples). In 
view of a possible mapping of the phase diagram of QCD 
using simulations of grand-canonical ensembles in W5 , it 
is of great interest to extend the analysis of [III [2(| about 
inhomogeneous condensates, to CP5. Moreover, in our 
opinion it is important to understand quantitatively how 
the numerical predictions of the model, namely Equa- 
tion (|23p . arc affected by varying the quark masses, and 



introducing a third flavor. We plan to report on these 
topics in the next future. 
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